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We use observations of the continuum radio emission at 1420, 820, 408, 34.5 and 22 MHz to esti- 
mate the mean brightness temperatures of the HB 21 supernova remnant (SNR) at five frequencies. 
We also presented mean spectral index of HB 21. The spectra of HB 21 are estimated for mean 
temperatures versus frequency for 1420, 820, 408, 34.5 and 22 MHz. We also presented T — T plots 
of three frequency pairs: between 1420-34.5, 1420-22, 34.5-22 MHz. We noticed flatter spectral 
indices at frequencies below 408 MHz. Probably this is due to the absorption by thermal plasma at 
low frequencies. 
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I. INTRODUCTION 

The radio emission from supernova remnants (SNRs) 
is generally understood to be synchrotron emission from 
relativistic electrons moving in magnetic fields. 

The HB 21 remnant is listed in Green's catalogue of 
Galactic supernova remnants ^i, i3] as G89.0-(-4.7. It has 
mixed morphology, e.g., shell-like in radio and center- 
filled in the X-ray As described in Kothes et al. [1], 
this is an old SNR evolving through the later stages of the 
Sedov phase, or early into the radiative (where a signifi- 
cant amount of the shock energy is being radiated away) . 
Results from Byun et al. [h] suggest that the unusual 
radio features and the central thermal X-ray enhance- 
ments of HB 21 might be the result of an interaction 
with molecular clouds. One form of evidences for the in- 
teraction between the SNR and the molecular cloud in 
HB 21 is the discovery of shocked molecular clumps by 
Koo et al. |6|. The center- filled, thermal X-ray emission 
is suggested to be caused by interaction with molecular 
clouds. 

Flux-calibrated images and spectra of the extended 
remnant HB 21 in the optical range, are presented in 
Mavromatakis, Xilouris & Boumis [7]. Filamentary and 
patchy structures were detected and appear to be corre- 
lated with the radio emission. 

Significant spatial variations in spectral index are 
found in HB 21 Q. He consider different physical mech- 
anisms for spectral index changes, including a detailed 
consideration of ionization losses in the dense molecular 
gas interacting with HB 21. 

One of the main goals of this paper is to provide ev- 
idence for Leahy's @ theoretical proposal. He suggests 



that if the spectra below 408 MHz show fiatter spectral 
indices then absorption by thermal plasma would be the 
preferred mechanism. We show the evidence of spectrum 
fiattening at the low radio frequencies in the supernova 
remnant HB 21. 



II. DATA AND METHOD 

Measured data at five frequencies are used: 1420 
Ik 820 [3, 408 [lii, 34.5 [13 and 22 MHz 
13{ . These surveys are available in electronic 
form in "Flexible Image Transport System" (FITS) 
data format, at the site of " Max-Planck- Institut 
fill Radioastronomie" (MPIfR) near Bonn, Germany: 

Using 



http : //www . mpif r-bonn . mpg . de/survey . html 
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this onhne Survey Sampler allows users to select a region 
of the sky and obtain images and data at different fre- 
quencies. The 1420-MHz Stockert survey H h as resolu- 
tion 0°.59, the 820-MHz Dwingeloo survey l°-2, the 
408-MHz all-sky survey [U 0°.85, the 34.5-MHz Gau- 
ribidanur survey [l^ 0°.7 and 22 MHz p"] 1°.7. The cor- 
responding observations are given at the following rates 
(measured data) for both I and b: at 1420 MHz, ^ at 

820 MHz, ^ at 408 MHz , ^ at 34.5 MHz and ^ at 22 
MHz. The effective sensitivities are about 50 mK T;, (Tf, 
is for an average brightness temperature), 0.20 K, 1.0 K, 
about 700 K and about 500 K, respectively. 

We extracted observed brightness temperatures from 
this data and then using programs in C and FORTRAN, 
obtained results shown in this paper. 

The area of HB 21 SNR is enclosed with brightness 
temperature contours. The maps of a region in HB 21 
SNR, in new Galactic coordinates (/, 6), with contours of 
the brightness temperatures Tb are plotted in Figs. [T][5] 
The contour lines correspond to the brightness temper- 
atures: minimum, maximum and nine contours in be- 
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TABLE I: Temperatures and brightnesses of HB 21 at 1420, 820, 408, 34.5 and 22 MHz. 



Frequency 


Temperature limits 


Temperature 


Brightness 


(MHz) 




(K) 


(10-^^ W/(m2 Hz Sr)) 


1420 


6.0, 9.0 


1.50 ± 0.05 


9.26 ± 0.30 


820 


16.5, 26.0 


4.91 ± 0.20 


10.13 ± 0.40 


408 


87, 140 


32.2 ± 1.0 


16.43 ± 0.50 


34.5 


49000, 81000 


18237 ± 700 


66.60 ± 2.56 


22 


105000, 144000 


28530 ± 500 


43.34 ± 0.76 



tween. The space between minimum and maximum con- 
tom'S defines HB 21 SNR. The Galactic longitude and 
latitude intervals for HB 21 SNR are the following: I — 
[90°, 87.5°], b = [3.5°, 6°]. We used the same method 
of calculation as given in Borka [3] for Galactic radio 
loops I- VI, Borka, Milogradov- Turin & Urosevic jlSll for 
Loops V and VI, Borka Jovanovic & Urosevic [lg| for 
Monoceros loop and in Borka Jovanovic & Urosevic [l3| 
for Gygnus loop. 

The mean temperatures and surface brightnesses of 
this SNR are computed using data taken from radio- 
continuum surveys at 1420, 820, 408, 34.5 and 22 MHz. 
The areas over which an average brightness temperature 
is determined at each of the five frequencies are taken to 
be as similar as possible within the limits of measurement 
accuracy. However, some differences between these areas 
still remain and we think that the major causes of dif- 
fering borders between the five frequencies are measured 
accuracy and small random and systematic errors in the 
data. The surface brightness of SNR must be above the 
sensitivity limit of the observations and must be clearly 
distinguishable from the Galactic background emission 
p^ . In order to evaluate brightness temperatures over 
the SNR we took into account background radiation (see 
Webster [l^). Borders enclosing the SNR are defined to 
separate the SNR from the background. 

We have subtracted the background radiation, in or- 
der to derive the mean brightness temperature of the 
SNR alone. The areas over which an average brightness 
temperature is determined at each of the five frequencies 
are taken to be as similar as possible within the limits 
of measurement accuracy. Tmin from Table [J means the 
lower temperature limit between the background and the 
SNR, and Tmax means the upper temperature of the SNR. 
We used all measured values inside the corresponding re- 
gions of I and 6, to calculate the brightness temperature 
of a SNR including the background. The mean bright- 
ness temperature for the SNR is found by subtracting the 
mean value of background brightness temperature from 
the mean value of the brightness temperature over the 
area of the SNR. 

After deriving the mean brightness temperatures T^, 
we derived surface brightnesses by: 

= (2fct/Vc')T,. (1) 

where k is Boltzmann constant and c the speed of light. 
Results are given in Table [J 



TABLE II: Brightnesses of HB 21 reduced to 1000 MHz, using 
spectral index a = 0.45 ± 0.07. 



Frequency 
(MHz) 


Brightness reduced at 1000 MHz 
(10^^^ W/(m2 Hz Sr)) 


1420 


10.83 ± 0.64 


820 


9.27 ± 0.24 


408 


11.01 ± 0.38 


34.5 


14.80 ± 3.10 


22 


7.88 ± 2.08 



III. DISCUSSION AND RESULTS 

The results for temperature limits are given in Table ID 
Tmin, given in the second column of Table HI is the lower 
temperature limit, while Tmax is the upper temperature 
limit given in the third column of Table ID There are 
some other sources near the SNR and they do not affect 
the calculation. These temperature limits enable us to 
distinguish the SNR from background and also from ex- 
ternal sources. Then we derived the temperatures and 
the surface brightnesses using equation (1) for each fre- 
quency. These results are given in the forth and fifth 
columns of Table HI respectively. Brightnesses of HB 21 
reduced to 1000 MHz, using spectral index a = 0.45 ± 
0.07 are given in Table HIl 

The radio map of HB 21 SNR at 1420 MHz with con- 
tours of brightness temperature is given in Fig. 1. This 
SNR has position: I = [90.°, 87.5°]; b = [3.5°, 6°], in 
new Galactic coordinates (l, b). Eleven contours plotted 
represent the temperatures Tmin and Tmax from Table H] 
and nine contours in between. The contour interval is 0.3 
K Tb, starting from the lowest temperature of 6 K up to 
9 K. Also, the corresponding temperature scale is given. 
In Fig. 2 is presented the radio map of HB 21 SNR at 
820 MHz with contours of brightness temperature. The 
contour interval is 9.5 K Tb, starting from the lowest tem- 
perature of 16.5 K up to 26 K. The radio map of HB 21 
SNR at 408 MHz with contours of brightness tempera- 
ture is given in Fig. 3. The contour interval is 5.3 K Tb, 
starting from the lowest temperature of 87 K up to 140 
K. In Fig. 4 is presented the radio map of HB 21 SNR at 
34.5 MHz with contours of brightness temperature. The 
contour interval is 3200 K Tb, starting from the lowest 
temperature of 49000 K up to 81000 K. The radio map of 
HB 21 SNR at 22 MHz with contours of brightness tem- 
perature is given in Fig. 5. The contour interval is 3900 
K Tb, starting from the lowest temperature of 105000 K 
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FIG. 1: Top: The radio map of HB 21 SNR at 1420 MHz with 
contours of brightness temperature. This SNR has position: 
I — [90.°, 87.5°]; b — [3.5°, 6°], in new Galactic coordinates 
{I, b). Eleven contours plotted represent the temperatures 
rmin and Tinax from Table |T] and nine contours in between. 
The contour interval is 0.3 K Tb, starting from the lowest 
temperature of 6 K up to 9 K. The corresponding temperature 
scale is given. Bottom: the 1420 MHz area map of Cygnus. 



up to 144000 K. The corresponding temperature scale is 
given for brightness temperatures of HB 21 at 1420, 820, 
408, 34.5 and 22 MHz at Figs. dB)-®, respectively. 
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FIG. 2: The same as Fig. 1, but for 820 MHz. The contour 
interval is 9.5 K Tb, starting from the lowest temperature of 
16.5 K up to 26 K. 



A. Spectrum 

The dominant emission mechanism from Galactic su- 
pernova remnants (SNRs) at radio frequencies is syn- 
chrotron emission. However in some SNRs, observations 
over a very broad range of radio frequencies reveal a cur- 
vature in the spectra of these sources @, [13] . The detec- 
tion of thermal emission at radio frequencies from Galac- 
tic SNRs shows us whether these sources are interacting 
with adjacent molecular clouds. It might also be useful 
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FIG. 3; The same as Fig. 1, but for 408 MHz. The contour 
interval is 5.3 K Tb, starting from the lowest temperature of 
87 K up to 140 K. 
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FIG. 4: The same as Fig. 1, but for 34.5 MHz. The contour 
interval is 3200 K Tb, starting from the lowest temperature 
of 49000 K up to 81000 K. 



for estimating the ambient density near SNRs using ra- 
dio continuum data (2l| . The presence of a significant 
amount of thermal emission will produce a curvature in 
the observed radio spectrum of the SNR, particularly for 
frequencies of 1 GHz and greater [21} . 

The radio continuum spectrum of HB 21 SNR (tem- 
perature vs frequency), is presented in Fig. |6l Spectra 
for five measurements - at 22, 34.5, 408, 820 and 1420 
MHz are given by solid line, spectrum obtained from the 
three frequencies - 408, 820 and 1420 MHz is given by 



dashed line, and spectrum obtained from the two lowest 
frequencies - 22 and 34.5 MHz is given by dash-dot line. 
The spectrum was generated using mean temperatures 
at five different frequencies (see Fig. |6]). The brightness 
temperature spectral index is defined by Equation 
All five frequencies 1420, 820, 408, 34.5 and 22 MHz, 
from linear fit, give — 0.45 ± 0.07. Frequencies 34.5 
and 22 MHz, lie on very low ends of the spectrum, as 
presented in Fig. Inland we analyze them also separately. 
Derived value a2 = -1.01 (from two frequencies), sug- 
gests possible thermal absorption at very low end of the 
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FIG. 5: The same as Fig. 1, but for 22 MHz. The contour 
interval is 3900 K Tb, starting from the lowest temperature 
of 105000 K up to 144000 K. 



spectrum. 

The curved spectrum of the thh'd order of HB 21, ob- 
tained from the five frequencies is presented in Fig. [T] 
Obtained spectra show that a curvature [22i,^] is present 
in the radio spectrum at low frequencies, see Figure [71 In 
the papers Baars, Mezger & Wendker [22], Milogradov- 
Turin & Nikolic ^23,] spectra are represented by second 
order polynomial fit. For our spectra we used polynomial 
fit of third order. This third order polynomial regression 
line fits properly temperature versus frequency at five 
measurements - at 22, 34.5, 408, 820 and 1420 MHz. 
The parameters of the third order polynomial fit are: a 
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FIG. 6: The radio continum spectrum of HB 21 SNR: tem- 
perature versus frequency. Spectrum for five measurements 
- at 22, 34.5, 408, 820 and 1420 MHz is given by solid line, 
spectrum obtained from the three frequencies - 408, 820 and 
1420 MHz is given by dashed line, and spectrum obtained 
from the two lowest frequencies - 22 and 34.5 MHz is given 
by dash-dot line. 




FIG. 7: The curved spectrum of the third order of HB 21, 
obtained from the five frequencies. 



= -0.582641 ± 0.1809 (31.04%), b = 9.38084 ± 0.2556 
(2.725%), c = -5.15063 ± 0.1133 (2.201%), d = 0.714124 
± 0.01609 (2.254%). 

We are using two equations: 



logT 



b log v 



c log^ + d log^ V. 



(2) 



(3) 



6 




200 400 600 800 1000 1200 1400 
V (MHz) 

FIG. 8; Variations of spectral index a''{i') with frequency in 
the case of HB 21. Spectral index a'^{i') is obtained from the 
third order polynomial fit presented in Fig. 7. Open circles 
represent values of spectral indices q'^(/^) at the five observed 
frequencies. 



Equation ^ represents the dependence between 
brightness temperature, frequency as weh as brightness 
temperature spectral index /3, (/3 = a + 2, where a rep- 
resents flux density spectral index) and if is a constant. 
Equation ([3]) represents third order polynomial fit and 
tells us about dependance T of v. Using equations ^ 
and ([3]) and after taking a log of both sides of equations 
and derivation of both sides with respect of we obtained 
expression: 



a^iv) = -{b + 2clogu + 3d log^ i^) 



(4) 



where index c denotes curved spectrum. 

We want to stress that the same procedure can be used 
if spectra is given in flux density because relation ([T]) 
holds between brightness and brightness temperature 
Ti, and between flux density Si, and brightness rela- 
tion holds: 



is obtained from the third order polynomial fit presented 
in Fig. [71 Open circles represent values of spectral in- 
dices at the five observed frequencies. We noticed fiat- 
ter spectral indices at low frequencies below about 200 
MHz. Probably this is due to the thermal mechanism 
at frequencies below about 200 MHz or below observed 
frequency of 408 MHz. Absorption by thermal plasma is 
probably the preferred mechanism in the low frequency 
range of the spectra as it is proposed by Leahy 

From the curved spectrum, we can also calculate spec- 
tral index a'^ using equation Q . In this case a is function 
of ly, and therefore for mean a we use: 



where i^min — 22 MHz and 
Using eqs. (Hj) and (O we obtain a' 



= 1420 MHz. 



(7) 



= 22 MHz and Vr, 
-0.97. 



= 0.46. When 
34.5 MHz the result 



22-1420 



we put Vmi' 
is <522_34.5 

Since we deal with mean brightness temperature, we 
lose the information about index variation because for 
each frequency we have only one average value for bright- 
ness temperature. Using equation ^ we get a connec- 
tion between brightness temperature and spectral index 
a. We plot the dependence logT vs log:^. The depen- 
dence log V between 408, 820 and 1420 MHz can be nicely 
represented by straight line. Leahy Q obtained mean 
spectral index a = 0.45 using observations presented at 
408 and 1420 MHz, and from our Fig. 8 using only points 
at 408 and 1420 MHz we obtained the a = 0.47 which is 
in excellent agreement with Leahy's value. But if we take 
into account low frequencies 22 and 34.5 MHz we can find 
that now dependence log T vs log v is nicely represented 
by the polynomial fit of third order. 

This indicates that the spectrum of HB 21 is a com- 
bination of synchrotron and a thermal components. The 
presence of the additional component in the radio spec- 
trum of HB 21 suggests that HB 21 is probably interact- 
ing with an adjacent molecular cloud. 



s.^^^n (5) 

where Q, is the solid angle. 

The new equation will be the same as equation ([3]), but 
with different coefficients. 

log Sy^a + y log v + c log^ v + d' log^ v. (6) 

We choose spectrum given in Tb because our original 
data are given in Tb, not in the flux density. 

Using these parameters and equation ^ we can derive 
the dependence of spectral index with frequency along 
the entire frequency range from 22 to 1420 MHz. Varia- 
tions of spectral index a'^iy) with frequency in the case 
of HB 21 is presented in Fig. [Sj Spectral index a^{v) 



B. T -T plot 

The measured data have different resolutions for differ- 
ent frequencies. In order to obtain T — T plots the data 
are retabulated. We convolved data at 1420, 34.5 and 22 
MHz to 0°.25 X 0°.25 resolution. Sampfing rates of the 
820 MHz and 408 MHz survey are much more crude and 
we did not take it into account. Then, for each frequency 
pair we used only the common points (with the same 
{I, b)) which belong to the loop area at both frequencies. 
In that way we reduced loop area to the same area for 
different frequencies. The obtained T — T plots for three 
pairs of frequencies enabled calculating the spectral in- 
dices. We calculated two a values for each of these tree 
frequency pairs: between 1420-34.5, 1420-22, 34.5-22 
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MHz. For each of the tree frequency pairs, by interchang- 
ing the dependent and independent variables we have ob- 
tained two a values for each pair and the mean value of 
these fit results is adopted as the radio spectral index, as 
suggested in Uyaniker et al. [l^l- Regarding two pairs 
of frequencies 1420-34.5, 1420-22 the average value of 
spectral index from T - T is < utt >i420-22, 1420-34. 5 
= 0.47 ± 0.29. Regarding the lowest pair of frequencies 
34.5-22, the average value of spectral index from T ~T 
is < UTT >34.5-22 = -0.82±0.82. 

It can be noticed that this value agrees well with the 
corresponding value obtained from spectrum, as expected 
(see Uyaniker et al. [13] )■ 

Leahy Q determined in HB 21 SNR a mean spectral 
index of 0.45. Our mean spectral index is in agreement 
with Leahy findings. He found significant spectral index 
variations and concluded that thermal absorption was 
the preferred mechanism fs'] . He concluded that if ther- 
mal plasma absorption or ionization losses are the cor- 
rect mechanism, at frequencies below 408 MHz should 
show flatter spectral indices and if one should be able 
to distinguish which mechanism is working, since ioniza- 
tion losses will have no indices flatter than 0.25 @. We 
obtained flatter spectral indices then 0.25 for the HB 21 
SNR in the low frequency domain, and we can expect that 
thermal plasma absorption is responsible for the spectral 
flattening at the lowest radio frequencies. 



about 500 K Tb for 22 MHz. 

We present the radio continuum spectrum of the HB 
21 SNR using average brightness temperatures at five fre- 
quencies. As it can be seen from Fig. [51 given linear fit 
provides reliable spectral index. Our analysis indicates 
that significant spectral variations for different frequen- 
cies in spectral index are found in HB 21. We noticed 
flatter spectral indices at frequencies below 408 MHz as 
it is proposed by Leahy Q. 

This indicates that the spectrum of HB 21 is a com- 
bination of synchrotron and thermal components. The 
presence of the additional component in the radio spec- 
trum of HB 21 suggests that this SNR is interacting with 
an adjacent molecular cloud. 

In Leahy Q, observations of the SNR HB 21 are pre- 
sented at 408 and 1420 MHz. He considered different 
physical mechanisms for spectral index variations. He 
concluded that if the spectra below 408 MHz show flat- 
ter spectral indices then absorption by thermal plasma 
is the preferred mechanism. Our obtained spectra at 5 
frequencies leads to a similar conclusion. 

We note an obvious flattening of the spectral indices. 
In future work one might investigate a broader frequency 
range. Probably, thermal emission at higher frequencies 
will produce a curvature in the radio spectrum of the HB 
21 SNR, particularly for frequencies higher then 1 GHz 




IV. CONCLUSIONS 

We use observations of the continuum radio emission at 
1420, 820, 408, 34.5 and 22 MHz for estimations of the 
mean brightness temperatures and surface brightnesses 
of the HB 21 SNR. The sensitivity of the brightness tem- 
peratures are: 50 mK for 1420 MHz, 0.2 K for 820 MHz, 
1.0 K for 408 MHz, about 700 K Tb for 34.5 MHz and 
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